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1.  INTRODUCTION 

Shock-boundary  layer  interaction  can  significantly  influence  not  only  the  local 
transonic  flow  on  missiles,  wings  and  turbine  blades  but  its  influence  can  also  ex¬ 
tend  downstream  within  the  boundary  layer  and  thereby  alter  the  global  aerodynamic 
properties  of  lift,  drag  and  pitching  moment.  It  is  therefore  important  that  shock¬ 
boundary  layer  Interactions  and  their  Reynolds  and  Mach  number-scaling  be  properly 
modeled  in  engineering  flow  field  prediction  methods  for  supercritical  aerodynamic 
bodies.  This  paper  describes  the  application  of  a  non-asymptotlc  triple  deck  theory 
of  transonic  shock-turbulent  boundary  layer  interaction  which  provides  such  a  tool 
for  non-separating  two-dimensional  flows  over  a  wide  range  of  practical  Reynolds 
numbers.  Section  2  contains  a  brief  description  of  the  essential  features  of  the 
theoretical  model.  Sections  3  and  4  then  describe  how  the  results  of  a  comprehen¬ 
sive  parametric  study  of  this  theory  may  be  used  to  develop  a  generalized  "viscous 
wedge"  model  of  the  local  interaction  which  embodies  proper  scaling  behavior  as 
well  as  an  approximate  account  of  incipient  separation  that  is  in  good  agreement 
with  experimental  trends.  In  Section  5  we  examine  the  application  of  this  theory 
as  an  element  in  global  viscous  flow  field  analyses  of  supercritical  airfoils.  In 
such  problems  it  will  be  shown  that  even  in  non-separating  cases  the  changes  across 
the  interaction  may  significantly  alter  the  subsequent  turbulent  boundary  layer  be¬ 
havior  for  appreciable  distances,  especially  when  larger  downstream  adverse  prersure 
gradients  are  present. 

2.  BRIEF  OUTLINE  OF  THE  LOCAL  INTERACTION  THEORY 

Unlike  significantly-separated  flow  where  the  disturbance  flow  pattern  as -or I 
with  a  nearly-normal  shock-boundary  layer  interaction  is  a  very  complicated  an^  in¬ 
volving  a  bifurcated  shock  pattern',  the  unseparated  case  pertaining  to  turbulent 
boundary  layers  up  to  roughly  Mj  «  1.3  has  instead  a  much  simpler  type  of  lp»era<.  tiun 
pattern  which  is  more  amenable  to  analytical  treatment.  With  some  judic <ous  simpli¬ 
fications,  it  is  possible  to  construct  a  fundamentally-based  approximate  theory  .  f 
the  problem  in  the  latter  case.  Consider  a  known  adiabatic  boundary  layet  pro  11  ><s 
Mgfy)  subjected  to  small  transonic  disturbances  due  to  an  impinging  weak  and  nearly 
normal  shock.  In  the  practical  Reynolds  number  range  of  interest  here  (105<Re^<  r), 
it  has  been  established2-4  that  the  local  interaction  disturbance  field  in  the  ne  gh  - 
borhood  of  the  impinging  shock  organizes  itself  into  three  basic  layered  regions  >-r 
"decks"  (Fig.  1):  1)  an  outer  region  of  potential  invlscld  flow  above  the  boundary 
layer,  rrtiich  contains  the  incident  shock  and  interactive  wave  systems;  2)  an  late' 
mediate  deck  of  frozen  shear  stress-rotational  lnviscid  disturbance  flow  occupy ;r.* 
the  outer  90X  or  more  of  the  incoming  boundary  layer  thickness;  3)  an  inner  shear 
disturbance  sublayer  adjacent  to  the  wall  which  accounts  for  the  interactive  skit 
friction  perturbations  (and  hence  any  possible  incipient  separation)  plus  most  of 
the  upstream  Influence  of  the  interaction.  The  "forcing  function"  of  the  problem 
here  is  thus  impressed  by  the  outer  deck  upon  the  boundary  layer;  the  middle 
couples  this  to  the  response  of  the  inner  deck  but  in  so  doing  can  Itself  modify 
the  disturbance  field  to  some  extent,  while  the  slow  viscous  flow  in  the  thin  Inrn 
deck  reacts  very  strongly  to  the  pressure  gradient  disturbances  Imposed  by  these 
overlying  decks.  In  treating  this  interactive  field  we  employ  a  non-asymptotlc  an*’. fod 
that  is  an  extension  to  turbulent  flow  of  Llghthill's  approach6  because  of  Its  es¬ 
sential  soundness  and  adaptability  to  practical  engineering  problems,  similarity  tc 
related  types  of  multiple-deck  approaches  that  have  proven  highly  successful  In  treat¬ 
ing  turbulent  boundary  layer  response  to  strong  adverse  pressure  gradients,  and  the 
large  body  of  turbulent  boundary  layer  Interaction  data  plus  recent  numerical  studic- 
wlth  the  full  Navler-Stokes  equations  which  support  the  predicted  results  in  a  varte*v 
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of  problems  (see  the  survey  in  Ref.  5).  Moreover,  this  approach  provides  at  realistic 
Reynolds  numbers  a  treatment  of  the  inner  deck  pressure  gradient  terms  plus  the 
middle  deck  3p/3y  and  streamline  divergence  effects,  along  vith  simplifying  ap¬ 
proximations  that  render  the  resulting  theory  tractable  from  an  engineering  stand¬ 
point  . 

A  very  detailed  description  of  the  above-mentioned  non-asymptotlc  triple  deck 
analysis  can  be  found  in  Refs.  5  and  7  and  hence  will  not  be  given  here.  The  re¬ 
sulting  predictions,  such  as  typically  illustrated  in  Fig.  2,  describe  all  the  es¬ 
sential  global  features  of  the  mixed  transonic  character  of  the  non-separating  norma) 
shock-turbulent  boundary  layer  interaction  problem  including  the  interactive  pressure 
distribution  and  upstream  influence,  displacement  thickness  and  local  shape  factor, 
and  interactive  skin  friction  up  to  incipient  separation.  This  interaction  theory 
employs  for  the  Incoming  turbulent  boundary  layer  velocity  profile  a  very  general 
Composite  Law  of  the  Hall-Law  of  the  Wake  profile  model  due  to  Walz8 ,  which  is 
characterized  by  three  parameters  (Mi,  boundary  layer  thickness  Reynolds  number  and 
the  incoming  shape  factor).  The  influence  of  both  shock  obliquity^  and  wall  curva¬ 
ture10  have  also  been  examined  in  detail  and  incorporated  into  the  theory.  Very 
extensive  parametric  studies  and  detailed  comparisons  with  experiment  have  shown 
that  it  gives  a  very  good  account  of  all  the  important  engineering  features  of  the 
interaction  over  a  wide  range  of  Mach  and  Reynolds  number  conditions.  Moreover,  the 
Important  but  heretofore-ignored  influence  of  incoming  boundary  layer  shape  factor 
Hj ^  (and  hence  the  upstream  pressure  gradient  history)  has  been  determined^*7. 

3.  A  GENERALIZED  VISCOUS-RAMP  MODEL  OF  THE  INTERACTION 

In  certain  engineering  applications  to  global  flow  field  analysis  computer  o»> 
grams  for  wings  or  turbine  blades,  it  has  proved  expedient  to  model  the  interaction 
as  a  simple  local  lnvlscld  6*  -  "bump"  or  "ramp."  A  serious  deficiency  of  this  ap¬ 
proach  is  that  it  does  not  account  for  the  dependence  of  the  bump  shape  aed  size  on 
Reynolds  number,  shock  strength  and  boundary  layer  shape  factor,  while  the  additional 
interaction  effects  on  the  downstream  boundary  layer  (such  as  Cf  reduction;  ire  ig¬ 
nored  altogether.  With  the  aforementioned  parametric  study  results  in  hand,  hovpver, 
the  present  theory  provides  a  much  improved 'Viscous  ramp"  representation  of  che  in 
teraction  (see  Fig.  3A)  whose  key  physical  features  have  the  correct  dependence  ci 
Mi ,  Re^  and 

Results  for  these  viscous  wedge  properties  taken  from  Ref.  11  are  presented  ;n 
Figs.  3B-3E,  where  the  upstream  and  downstream  influence  distances,  the  slope  and 
overall  height  of  the  6*  -  bump  are  plotted  along  with  the  downstream  Ct/CfQ  values 
that  may  be  needed  along  with  6*  to  re- initialize  a  subsequent  turbulent  boundary 
layer  calculation  downstream.  We  note  in  general  that  the  overall  scale  of  the  to 
teraction  (which  can  be  a  sensitive  effect  in  both  steady  and  unsteady  applications 
where  such  viscous  wedge  models  are  employed)  does  not  scale  according  to  the  un¬ 
disturbed  boundary  layer  thickness  even  in  the  non-separating  case.  It  is  further 
noted  that  the  viscous  wedge  slopes  are  in  rough  agreement  with  the  maximum  attached 
shock  deflection  value  observed  empirically12,  although  here  of  course  there  Is  a 
dependence  on  Re^  and  Hj.  as  well  as  Mach  number.  Finally,  in  all  of  these  curves 
we  see  a  significant  dependence  on  the  incoming  boundary  layer  shape  factor  that  wouJd 
appear  to  be  an  Important  consideration  in  practical  applications. 

Closed-form  analytical  fits  to  the  various  curves  of  Fig.  3  have  been  deveJope<>  ‘ ' 
which  provide  a  very  rapid  yet  complete  modeling  of  the  local  Interaction  effects  to- 
incorporation  as  a  local  module  in  global  lnviscid-boundary  layer  analysis  or  design 
programs  for  supercritical  airfoils. 
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4.  INCIPIENT  SEPARATION 

The  present  theory,  although  it  breaks  down  at  separation,  does  yield  a  useful 
indication  of  incipient  separation  where  Cf  .  0,  owing  to  the  particular  attention 

paid  to  the  treatment  of  the  local  interactive  skin  friction  behavior  •  Since  this 
indication  is  of  great  practical  interest,  a  parametric  study  of  incipient  separa¬ 
tion  conditions  inherent  in  the  present  theory  was  carried  out;  the  results  for  a 
normal  shock  on  a  flat  surface  are  presented  in  Fig.  4A  where  the  shock  Mach  number 
above  which  incipient  separation  occurs  is  plotted  as  a  function  of  the  Reynolds  num¬ 
ber  with  the  shape  factor  as  a  parameter;  also  shown  in  the  Figure  is  the  approximate 
experimental  boundary  determined  by  a  careful  examination11  of  a  large  number  of  trans¬ 
onic  Interaction  tests,  besides  Nussdorfer*  s1  ®  M  ~1.30  criterion  for  turbulent  flow.  It 
is  seen  that  the  theoretical  prediction  of  a  gradual  increase  in  the  incipient  separa¬ 
tion  Mach  number  value  with  Reynolds  number  is  in  agreement  with  the  trend  of  the 
data;  moreover,  the  theoretical  prediction  of  only  a  small  influence  of  shape  factor 
on  the  incipient  separation  conditions  is  also  borne  out  by  the  lack  of  any  consistent 
H-ef feet  for  the  name  Re  discernible  in  the  data  (Squire  has  observed  a  similar  insen¬ 
sitivity  to  Hj.  in  purely  supersonic  flow  interactions14).  The  absolute  values  of 
the  incipient  separation  Mach  number  predicted  by  the  present  interaction  theory  are 
seen  to  be  consistently  slightly  lower  than  the  average  experimental  value;  this  is 
attributable  to  the  combined  effects  of  the  linearized  inner  deck  theory  (which  over- 
predicts  the  pressure  gradient  effect  on  Cf)  and  the  assumption  of  a  normal  shock 
when  in  fact  most  of  the  experiments  likely  entail  some  shock  obliquity  (which  also 
delays  separation  to  somewhat  higher  shock-strengths).  It  is  interesting  to  note 
that  Nussdorfer' s1 3  original  incipient  separation  criterion,  based  as  it  was  on  a 
very  limited  base,  does  roughly  go  through  the  average  of  the  data  although  it  does 
not  account  for  the  proper  Reynolds  number  effect. 

Fig.  4B  shows  the  influence  of  wail  curvature;  it  is  seen  to  have  onlv  a  small 
effect  in  delaying  incipient  separation  to  a  slightly  higher  shock  strengtr>  for  » 
given  Re^  and  Hf,  this  being  of  the  same  order  as  the  experimental  data  b„r.d . 

5.  APPLICATION  TO  GLOBAL  FLOW  FIELD  ANALYSIS 
5.1  Downstream  Effects  from  a  Shock-Boundary  Layer  Interaction 

In  addition  to  the  Increased  displacement  thickness,  the  foregoing  discussion 
shows  that  the  skin  friction  level  following  the  interaction  is  significantly  redu.-d; 
combined  with  the  attendant  sensitivity  to  the  profile  shape,  this  suggests  that 
subsequent  downstream  boundary  layer  development  may  retain  a  memory  of  the  Intet ac¬ 
tion  effects  for  a  considerable  distance  (over  and  above  a  simple  thickening),  par 
ticularly  as  regards  possible  incipient  separation  in  any  adverse  pressure  giadleo- s 
downstream.  As  indicated  in  Fig.  5  this  downstream  "interaction  after-effect"  in 
the  boundary  layer  Influences  the  sensitive  trailing  edge  region  and  thus  may  be  ia- 
portent  in  the  design  and  analysis  of  rear-loaded  airfoils,  especially  at  higher  ltn 
coefficients  with  increasingly-aft  shock  locations;  it  likewise  may  be  Important  on 
three  dimensional  wing  configurations  where  the  shock  interaction  zones  are  well  «fr. 

The  aforementioned  after-effect  question  was  therefore  subjected  to  detailed 
study  using  the  two-layer  turbulent  boundary  layer  program  of  Mosesls  as  a  model  ct 
the  downstream  viscous  flow;  the  program  is  coupled  to  the  present  interaction  theot> 
by  initializing  it  with  the  post-interactive  flow  properties  so  as  to  account  fully 
(both  Cf  and  4*),  partially  (&*  only)  or  not  at  all  for  the  changes  across  the  inter¬ 
action.  Calculations  were  then  made  of  the  subsequent  downstream  turbulent  bounder* 

behavior  (8,  C-,  8*,  6*)  in  various  constant  post-shock  adverse  pressure  gradient* 
typical  of  airfoils  for  different  assumed  local  interactive  shock  strengths  and  post 
tions  or  Reynolds  numbers.  The  results  serve  as  a  paradigm  of  the  downstream  sensi¬ 
tivity  question  in  real  flows. 

A  variety  of  cases  were  studied16,  typical  results  of  which  are  presented  in 


Fig.  6  where  we  show  Che  predicted  behavior  of  the  boundary  layer  shape  factor  and 
skin  friction  in  three  lncreasingly-strong  adverse  pressure  gradients  downstream 
of  an  interaction  occurring  at  a  typically  rearward  position;  the  consequences  of 
fully,  partially  or  negligently-treating  the  boundary  layer  changes  across  the  inter¬ 
action  are  indicated.  Generally,  it  is  seen  that  the  downstream  behavior  of  the 
boundary  layer  is  indeed  sensitive  to  detailed  modeling  of  the  interactive  effects 
and  that  this  sensitivity  increases  with  the  strength  of  the  downstream  adverse 
pressure  gradient.  The  adverse  pressure  gradient  magnifies  the  subsequent  influ¬ 
ence  of  the  skin  friction  (as  well  as  the  £*-rise)  across  the  Interaction  so  that 
downstream  separation  tends  to  occur  earlier  than  would  be  predicted  by  either  neg¬ 
lecting  or  treating  only  the  6*  effect  of  the  upstream  interaction.  As  shown  in 
Fig.  7  ,  these  predictions  are  supported  by  a  comparison  with  boundary  layer  measure¬ 
ments  downstream  of  a  non-separating  shock  interaction  zone  on  a  supercritical  airfoil; 
both  the  skin  friction  and  shape  factor  data  are  poorly  predicted  when  the  interaction 
is  neglected  but  are  reasonably  well  predicted  when  the  complete  interaction  effects 
are  taken  into  account. 

Examination  of  many  such  results  leads  to  the  further  conclusion  that  such  in¬ 
teractive  after-effects  extend  at  least  20-302  chord  distances  downstream  on  a  typical 
airfoil  or  wing  and  increase  (as  expected)  with  either  larger  shock  strength  or 
decreasing  Reynolds  number.  If  the  trailing  edge  region  lies  within  this  range  of 
the  shock,  it  is  thus  seen  that  a  simple  thickening  effect  alone  is  not  sufficient 
to  account  for  the  interaction  and  may  result  in  an  inaccurate  prediction  of  the 
rearward  boundary  layer  shape  factor,  skin  friction  and  incipient  separation  proper¬ 
ties  including  their  scaling.  This  is  of  practical  importance  for  two  major  reasons: 

1)  in  regions  of  sustained  adverse  pressure  gradient  that  often  follow  the  short- 
scale  interaction  zone,  the  shape  of  the  velocity  profile  and  streamwise  shear  stress 
distribution  (as  well  as  thickness)  are  of  considerable  importance  to  the  aerodynamic 
design  of  an  airfoil  or  wing;  2)  the  altered  boundary  layer  properties  (especially 
possible  incipient  separation)  near  the  trailing  edge  and  into  the  wake  can  further 
exert  a  powerful  effect  on  the  overall  aerodynamics  via  their  influence  of  the  Kutta 
condition17  and  on  possible  buffet  onset. 

5.2  Supercritical  Wing  Section  Flow  Fields 

Nandanan  et  al12  have  carried  out  an  even  more  detailed  study  of  interactions  on 
actual  supercritical  airfoils  including  experimental  comparisons.  They  developed  a 
global  computational  method  for  transonic  airfoil  flow  analysis  which  incorporates  the 
present  analytical  solution  for  near-normal  shock-boundary  layer  Interaction  into  a 
state-of-the-art  viscous- lnvisc id  computation  code.  Theoretical  results  obtalneo 
with  this  method  were  compared  to  representative  data  from  boundary  layer  and  surface 
pressure  measurements  on  three  transonic  airfoils  in  the  DFVLR-AVA  (Gottingen)  Trans¬ 
onic  Hind  Tunnel;  some  examples  of  these  comparisons  are  shown  in  Figs.  8A  and  8B  The 
agreement  between  theory  and  experiment  in  both  the  boundary  layer  displacement  chick 
ness  and  the  surface  pressure  distributions  was,  for  all  test  cases  considered,  quite 
good.  The  associated  predictions  of  the  local  skin  friction  variation  through  the 
interaction  zone  also  agree  reasonably  well  with  the  values  Inferred  from  the  experi¬ 
mental  boundary  layer  profiles  via  the  Ludwig-Tlllman  relation  (see,  e.g..  Fig.  9  ) . 

The  results  of  this  investigation  indicated  that  treating  the  shock-boundary  layer 
interaction  by  conventional  boundary  layer  theory  generally  leads  to  a  slight  underpve 
diction  of  the  displacement  thickness  immediately  downstream  of  the  shock  and,  due  ra 
the  amplifying  effect  of  the  sustained  rear  adverse  pressure  gradients,  to  an  appre¬ 
ciable  underestimation  of  the  displacement  thickness  at  the  trailing  edge  (see  Fig. 10). 
The  latter  is  also  clearly  reflected  in  the  pressure  distributions  and  aerodynamic 
coefficients  compared.  Considering  these  results,  one  may  conclude  that  it  is  general  y 
necessary  to  include  a  physically  correct  treatment  of  shock  wave-boundary  layer  intei 
action  in  the  analysis  of  transonic  airfoil  flow. 
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6.  CONCLUDING  REMARKS 


The  results  of  this  study  have  shown  that  It  Is  now  possible  to  Incorporate 
as  an  Interactive  module  within  a  global  flow  field  analysis  the  correctly-modeled 
(and  scaled)  local  shock-boundary  layer  interaction  effects  for  the  non-separating 
case.  The  non-asymptotic  triple  deck  interaction  theory  involved  covers  a  wide 
range  of  practical  Reynolds  numbers  and  turbulent  boundary  layer  profile  shape  fac¬ 
tors  including  the  effect  of  wall  curvature;  moreover,  it  gives  an  approximate  in¬ 
dication  of  when  incipient  separation  occurs.  It  was  further  shown  that  such  theory 
is  generally  desirable  when  accurate  predictions  are  desired  in  the  important  trail¬ 
ing  edge  region  of  rear-loaded  supercritical  airfoils  because  the  detailed  changes 
across  an  upstream  interaction  can  significantly  alter  the  subsequent  turbulent 
boundary  layer  behavior  for  appreciable  distances  downstream. 
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